Theory of shot noise in single-walled metallic carbon nanotubes weakly coupled to 

nonmagnetic and ferromagnetic leads 
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We present theoretical study of shot noise in single wall metallic carbon nanotubes weakly cou- 
pled to either nonmagnetic or ferromagnetic leads. Using the real-time diagrammatic technique, we 
calculate the current, Fano factor and tunnel magnetoresistance in the sequential tunneling regime. 
It is shown that the differential conductance displays characteristic four-fold periodicity, indicating 
single-electron charging. Such a periodicity is also visible in tunnel magnetoresistance of the sys- 
tem as well as in the Fano factor. The present studies elucidate the impact of ferromagnetic (vs. 
nonmagnetic) contacts on the transport characteristics under consideration. 

PACS numbers: 72.25.Mk, 73.63.Kv, 85.75.-d, 73.23.Hk 
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I. INTRODUCTION 

In connection with more and more demanding minia- 
turization requirements of the modern electronics, there 
has recently been much interest in electric transport 
in low-dimensional nanostructured materials hke quan- 
tum dots and quantum wires (both conventional^^ and 
magnetio^i^i^'i'^) . Similarly, carbon nanotubes (CNTs) 
have been attracting wide attention for more than one 
and a half decade due to their exceptional transport 
and mechanical properties^ In particular, depending on 
how a graphite monolayer (graphene) is rolled up into 
a cylinder, the resulting nanotube may be either metal- 
lic or semiconducting. CNTs may be used as intercon- 
nects as well as switching and nonvolatile memory de- 
vices in nanoelectronicsjiSiii Furthermore, there is no 
doubt that CNTs are also quite promising for spintronic 
applications, because when placed between ferromag- 
netic electrodes they show a considerable giant (tunnel) 
magnetoresistance, GMR (TMR), ranging from a few, 
up to several tens percent ^'i^^'i^^iii^'ia^^^^^ 
Moreover, recent experiments on single wall CNTs with 
moderate coupling to external leadsi^ have shown that 
TMR in nanotubes can become negative. This behav- 
ior was accounted for theoretically in terms of spin- 
dependent phase shifts at interfacesi^ Therefore our the- 
oretical analysis of CNTs is extended so as to include, on 
the one hand detailed intrinsic energy-level structure of 
the given CNT (singlet vs. triplet energy level occupan- 
cies, etc.), and on the other hand the effect of magnetic 
electrodes on transport properties. It is now well-known 
that at low temperature the Coulomb blockade regime 
may manifest itself even when the so-called addition en- 
ergy's less than energy-level spacings. As shown in Refs. 
llHl], the conductance spectra of CNTs in that case 



reveal a characteristic four-fold structure with relatively 
high peaks (of the order of /h). 

In this paper we analyze theoretically the charge and 
spin transport in single wall CNTs weakly coupled to 
metallic leads - either nonmagnetic or ferromagnetic. We 



consider Coulomb blockade transport regime, where the 
mean-field model of Oreg et al3i is applicable. That is 
the case when the charging energy-to-mean level spac- 
ing ratio is small enough (and consequently the Lut- 
tinger parameter g in not drastically less than 1) - 
otherwise, the Luttinger-liquid model would be more 



appropriate 
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Our considerations are based on the 



real-time diagrammatic technique. Assuming realistic 
parameters of the system^^i^ we calculate the basic 
transport characteristics, i.e. the current, conductance, 
and TMR in the case of ferromagnetic contacts. Apart 
from this, we also calculate the shot noise (associated 
with discrete nature of charge) and show that the corre- 
sponding Fano factor F is suppressed below the Schot- 
tky value, F <1, beyond the Coulomb blockade regime. 
More specifically, in the transport regions where sequen- 
tial contribution to the current is dominant, the shot 
noise is sub-Poissonian with F slightly above 0.5. Since 
our considerations are limited to the first-order trans- 
port with respect to the coupling parameter, the above 
conclusion does not apply to regions where sequential 
tunneling processes are exponentially suppressed. It 
is noteworthy that the linear and nonlinear transport 
across carbon nanotubes has been recently addressed by 
Mayrhofer and Grifonii '^^i^^ These considerations how- 
ever deal with the case of (i) nonmagnetic leads, and 
(ii) basically large diameter nanotubes, for which both 
the exchange effects and back scattering processes can 
be neglected. Spin-polarized transport for coUincar and 
noncoUinear alignments of the leads' magnetic moments 
was discussed in Refs. [23|33j. However, the consider- 
ations presented in Ref. [33[ are again subjected to the 
constraints (ii) , whereas Ref. [1^ concerns only the linear 
response regime (zero bias limit). In the present paper 
we focus mainly on the shot noise, using a nonequilibrium 
approach applicable, in principle, to carbon nanotubes of 
arbitrary diameter. It seems that, at present, theoreti- 
cal studies on the shot noise in CNTs lag behind recent 
experimental ones (e.g. [34.35,36,]). On the other hand, 
we are not aware of any investigations (either theoretical 



2 



or experimental) on the shot noise in ferromagnetically 
contacted CNTs. Therefore, the main objective of this 
study is to fill, at least partially, these gaps. In addi- 
tion, we also believe that our considerations will stim- 
ulate experimental research on CNTs weakly coupled to 
ferromagnetic leads, and will be helpful in understanding 
future experiments. 

The paper is organized as follows: In Sec. II the model 
and computational method are shortly described. Sec- 
tions III and IV present results on electronic transport 
(including the shot noise) in the sequential tunneling 
approximation for nonmagnetic and ferromagnetic con- 
tacts, respectively. Finally, Sec. V provides a short sum- 
mary and conclusions. 



II. MODEL AND METHOD 

The system considered in this paper consists of a sin- 
gle wall metallic CNT and two electron reservoirs (elec- 
trodes) which are weakly coupled to the CNT. The elec- 
trodes can be either nonmagnetic or ferromagnetic. In 
the latter case, the net spin moments (and magnetiza- 
tions) of the leads are assumed to be collinear, i.e., they 
can form either parallel or antiparallel magnetic config- 
uration. The Hamiltonian H of the system takes the 
general form 



H = Hl + Hr + Hcnt + Ht . 



(1) 



The first two terms describe noninteracting itinerant elec- 
trons in the leads. 



(2) 



kcr 



for the left (r = L) and right (r = R) lead, with e^kcr 
being the energy of an electron with the wave vector k 
and spin a in the lead r, and c]^^.^ (crka) denoting the 
respective creation (annihilation) operator. 

The third term, -ffcNT, in the Hamiltonian, Eq. ([T]), 
describes the CNT and is assumed in the form introduced 
by Oreg et ali^ 



H, 



CNT 



where 



= diiad-tija, and dl.^ {dfj.ja) is the creation 



(annihilation) operator of an electron with spin a on the 
jth level in the subband /i {fJ, = 1,2). The discrete struc- 
ture of the energy levels results from the quantization of 
electron motion along the CNT, which is an appropri- 
ate starting point in the limit of weak coupling between 
the nanotubc and the leads. The corresponding energy 
Efij of the jth discrete level in the subband fi is given by 



Efj^j = jA + {fj, — 1)5, where A is the mean level spac- 
ing and 6 describes the energy mismatch between the 
discrete levels corresponding to the two subbands. The 
second term in Eq. ^ stands for the electrostatic energy 
of a charged CNT, with U denoting the charging energy 
and A^o being the charge on the nanotube induced by 
gate voltages. The next term corresponds to the on-level 
Coulomb interaction with 6U being the relevant on-site 
Coulomb parameter. Finally, the last term in Eq. ([3]) 
describes the exchange energy, with J being the relevant 
exchange parameter. Furthermore, the system is sym- 
metrically biased and we assume equal capacitive cou- 
plings to the left and right leads, so the dependence of 
the nanotube energy levels on the bias voltage may be 
neglected. The exchange effects described by J play an 
important role for small diameter nanotubcs (up to 2 nm 
or so) , whereas for nanotubes of larger diameters these ef- 
fects become negligible. Thus in the latter case transport 
properties can be qualitatively described without taking 
into account the exchange interactioni^i^^ Incidentally 
we show in the following that J happens to be rather 
small in the present model, too. We also notice that the 
mean- field model of Oreg et ali^, Eq. (3), can be used 
to analyze the ground state properties and particle-hole 
excitations of single wall carbon nanotubes, however it 
fails to describe collective excitations such as spin and 
charge density waves, spin-charge separation, etc. (Lut- 
tinger liquid behavior) 

The last term of Hamiltonian, H^, takes into account 
tunneling processes between the CNT and electrodes, 
and is assumed in the form, 



(4) 



■r— L.R k ^j(7 



where denotes the tunnel matrix elements between the 
lead r and the jth level (assumed to be spin-independent 
also for ferromagnetic leads). The jth level coupling to 
external leads can be described by Tj.j = 2TT\trj\'^Pr, 
with p'^ being the density of states in the lead r for 
spin a. The role of ferromagnetic leads is taken into 
account just via the spin-dependent density of states 
p^. By introducing the spin polarization pr of the lead 



r, Pr 



- fn+ 



(p+ — p^ )/{pt + Pr )j the couphng parame- 
ters FJTj- can be expressed as F,^-^"'' = F,.j(l ±Pr), with 
Trj = (r+ + F" )/2. Here, F+ and F;^. describe the 
coupling of the jth level to the spin-majority and spin- 
minority electron bands of the rth lead, respectively. 
When the leads are nonmagnetic, then F^- 



F-. In 

the following we assume Trj = F/2 for all values of the 
indices j and r. 

We analyze electronic transport through the single wall 
metallic CNT described by the Hamiltonian ([3]) , and the 
considerations are limited to the sequential (first order 
in tunneling processes) tunneling regime. The first-order 
tunneling gives the dominant contribution to the charge 
and spin currents for voltages above a certain threshold 
voltage, and is exponentially suppressed in the Coulomb 
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blockade regime. The effects due to higher-order tunnel- 
ing processes, e.g. cotunneling^^i^i^S are not taken into 
account. When the bias voltage exceeds the threshold 
voltage, the current flows due to tunneling of electrons 
one by one through the discrete energy levels of the nan- 
otube. In order to make the analysis most realistic, in 
numerical calculations wc have taken into account up to 
six different orbital levels, three in each subband, which 
results in 4^ many-body nanotube states 

Transport is calculated with the aid of the real-time 
diagrammatic technique i^i^i^ This technique consists 
in a systematic perturbation expansion of the reduced 
density matrix and current operator in the couplings 
Trj. To determine the stationary occupation probabil- 
ities, charge current, and shot noise in the sequential 
tunneling regime, we employ the matrix approach de- 
veloped in Ref. [41|, and introduce the respective self- 
energy matrices: W, W', W''-. The matrix W con- 
tains self-energies with one arbitrary row xo replaced by 
(r, . . . , r), which is due to the normalization of the prob- 
abilities, J2x^x ~ ^- ^^'^ elements W^^i of the matrix 
W describe the first-order tunneling transitions between 
the Ix) and \x') many-body states. They are given by^^ 



^xx'+^xx'^ where 



Kx' = ^nY^pUMs^-ex^) 



for X 7^ X' and W^^^ = -J^x'^x^x'x^ "^^^^ ^^^^^ = 
ly/jgle-Mrj/feBT _|_ 2] and p,r being the electrochemical po- 
tential of the lead r. The second matrix, W', denotes 
the full self-energy matrix with one internal vertex, as- 
sociated with the expansion of the tunneling Hamilto- 
nian, replaced by the current operator. Finally, the 
third matrix, W'', consists of self-energies with two in- 
ternal vertices replaced by the current operator. The 
current operator / is defined as / = (/r — /l)/2, with 

ir = -i{e/h)Y,^^Y.^,j{irjcl^JMa~Kjdl.^Cr],a) be- 
ing the current flowing from the nanotube into 
the lead r. The elements of the matrices W' 
and W'-'- can be expressed in terms of W^x' as^ 

Kx' - [Q{Nx'-Nx)~mx~Nx')]{W^x--Kx') 
and W^-^, = (1 — 25xx')Wxx' / respectively, where 
Nx = ^^ja and @{x) is the step function. 

Having calculated the respective matrices, the station- 
ary probabilities can be determined from the following 
equation 



41 



{WV)^^T5xxo. 



(5) 



where P is the vector containing the occupation proba- 
bilities. In turn, the sequential current flowing through 




FIG. 1; Linear conductance as a function of the gate voltage. 
The parameters are: A = 8.4 meV, f//A = 0.26, J/A = 0.12, 
(5C//A = 0.04, 5/A = 0.27, feT/A = 0.025, pl = pr = 0, 
X = 0.14, and F = 0.2 meV. 



the CNT can be calculated from 



/=-Tr{W^P}, 



(6) 



with Tr{A} denoting the trace of the matrix A. On 
the other hand, the zero-frequency shot noise, S = 

2 j-oo {mm + mm - , is given by 



S = — Tr{W"P + W'PW'P} , (7) 

h 

where P is determined from the equation WP = Q, with 

Qx'x ~ i^x' ~ ^x'x) (1 ~ ^x'xo)'~ 

Having found the current / and the zero-frequency 
shot noise S, one can determine the Fano factor F, 
F = 5/(261/1). The Fano factor describes the deviation 
of S from the Poissonian shot noise given by Sp = 2e|/|. 
Now we proceed to numerical results based on the for- 
malism described briefly above, and we start with the 
case of nonmagnetic leads. 



III. CARBON NANOTUBES COUPLED TO 
NONMAGNETIC LEADS 

In this section we present numerical results on elec- 
tronic transport through CNTs contacted to nonmag- 
netic leads. The parameters assumed to model electronic 
structure of the CNT-based quantum dot weakly coupled 
to the leads are as follows: A ~ 8.4 meV, U/A ~ 0.26, 
J/A = 0.12, 5U/A = 0.04, 5/A = 0.27, and F = 0.2 
meV. These parameters are comparable to those derived 
from experimental observationsi ^^i^^ Apart from this, we 
assume k-oT/ A — 0.025, which also corresponds to the 
temperature of the above mentioned experiments. In 
addition, to fit our numerical data to the experimental 
ones presented in Refs. [HIl^, we assume equal capac- 
itive couplings of the CNT to the left and right leads. 



Cl = Cr, = C, while for the gate capacitance we assume, 
Cg = xC, with X = 0.14. 

By applying a gate voltage one can continuously shift 
position of the discrete energy levels up or down, which 
leads to peaks in the linear conductance each time the 
discrete level of the CNT crosses the Fermi level of the 
leads. This can be seen in Fig. [1] which displays the lin- 
ear conductance vs the gate voltage. First of all, one can 
see that the conductance spectrum exactly reproduces 
the four-peak periodicity observed experimentally in the 
weak coupling regimci^i^ For the parameters assumed 
here the dominant energy scale is set by the level sepa- 
ration A within individual electron subbands. The other 
important energy parameters are the charging energy U 
associated with an extra electron on the nanotube and 
the mismatch between the two subbands 5. The remain- 
ing parameters, i.e. the exchange coupling J and on-level 
Coulomb correlation 5U , are smaller and their influence 
on the spectra is less pronounced, although still remark- 
able. 

Each peak in the conductance spectrum corresponds 
to the addition of a new electron to the system, which 
fills in a subsequent empty energy level of the nanotube. 
Since there are two electron subbands in CNTs, shifted in 
energy by 5, and the nanotube levels are spin degenerate, 
four electrons can be added to the system before another 
discrete level in the same electron subband (separated by 
A) can be occupied, and a new four-peak sequence in the 
linear conductance can emerge. Conductance between 
the peaks is strongly suppressed, which is particularly 
well pronounced for the regions between successive four- 
peak patterns. The first peak in each four-peak sequence 
corresponds to the situation when a new energy level be- 
comes occupied by a single electron. The second peak, 
in turn, corresponds to the case when this level becomes 
occupied by another electron of opposite spin orienta- 
tion. The next two peaks appear when the energy level 
in the second subband (slightly higher in energy due to a 
nonzero mismatch parameter) becomes occupied first by 
a single electron and then by two electrons of opposite 
spins. Such scenario holds for the assumed parameters, 
i.e. for J + 5U < S. When, however, the exchange pa- 
rameter becomes larger, J -|- SU > S, the triplet state 
may be formed and a new scenario of the level filing may 
appear 

Differential conductance in both the linear and non- 
linear response regimes is shown in Fig. [21 which clearly 
reveals the blockade regions (the black diamonds) and 
discrete electronic structure. Figure[T]can be thus consid- 
ered as a cross-section of Fig.[2]along the line correspond- 
ing to V = 0. The three large diamonds in Fig. [5] corre- 
spond to the plateaus between the successive four-peak 
patterns from Fig.[TJ The smaller diamonds, in turn, cor- 
respond to the blockade regions between the peaks inside 
the corresponding four-peak sequence. 

The Fano factor, calculated as a function of the gate 
voltage and for the same parameters as in Fig. [1] is 
shown in Fig. [3] for several values of the bias voltage. 
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FIG. 2: (color online) Differential conductance as a function 
of bias and gate voltages. The parameters are the same as in 
Fig. [I 



\eV\ > k-QT. The latter condition ensures that the 
shot noise dominates over the thermal Nyquist- Johnson 
noise4^ When \eV\ < k^T, the contribution to the noise 
due to thermal fluctuations is significant, and in the 
strictly linear response regime, y — > 0, the correspond- 
ing Fano factor diverges due to nonzero thermal noise 
and vanishing of the average current density. Let us look 
first at the gate voltage dependence of the Fano factor for 
the lowest bias voltage shown in Fig. [31 = 2 mV. The 
Fano factor fluctuates then between values correspond- 
ing to the sub- and super-Poissonian noise. The super- 
Poissonian shot noise (F > 1) occurs in the blockade 
regions, whereas the sub-Poissonian shot noise (F < 1) 
is observed outside the blockade regions, where the se- 
quential transport is energetically allowed. The minima 
in Fig. [3l correspond then to positions of the peaks in 
Fig. m and vice versa the maxima in Fig. [3] correspond 
to the regions between the peaks in Fig. [H It is note- 
worthy that this type of correlation strongly resembles 
that reported for a zero-dimensional conventional (InAs) 
quantum dots 4^ 

When the bias voltage increases (see Fig. [3] for y = 4 
mV), the shot noise becomes sub-Poissonian, except for 
those gate voltage regions, where sequential transport is 
still suppressed. This takes place only for gate voltages 
around the centers of the large diamonds from Fig. [2] 
When the bias voltage increases further (see the curves 
for y = 8 mV and V = 12 mV in Fig. [3]), the system 
is beyond the blockade regime independently of the gate 
voltage. The corresponding Fano factor is always smaller 
than unity, although it fluctuates with increasing the gate 
voltage. The suppression of the Fano factor below F = 1 
in the sequential transport regime is a consequence of the 
Coulomb and exchange correlations in electronic trans- 
port. This behavior is different from that observed in 
the blockade regime, where F > 1.— However, one should 
bear in mind that the sequential transport is negligible in 
the blockade regions and the current is dominated there 
by higher-order (e.g. cotunneling) contributions. Thus, 
the Fano factor in the blockade regions may be signif- 
icantly modified (in particular reduced) by the cotun- 
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FIG. 3: Fano factor as a function of gate voltage for several 
values of bias voltage V. The parameters are the same as in 

Fig.m 



neling current. This, however, is beyond the scope of 
this paper, in which the considerations are limited to the 
first-order tunneling processes. 

The bias dependence of the Fano factor is explicitly 
shown in Fig. 2] for two different values of the gate volt- 
age. The first value (left column in Fig.Hl), Vg ~ —0.2666 



The first value (left column in Fig.d]), Vg 
V, corresponds to the center of the small diamond shown 
in Fig. [21 whereas the second value (right column in 
Fig. El), Vg ~ —0.1202 V, corresponds to the center of 
the large diamond. In order to see a clear correlation 
with transport properties, the current and differential 
conductance are also shown there. Firstly, the regions 
where the current is suppressed are clearly visible, while 
for voltages above the threshold voltage electrons tunnel 
one by one through the system giving rise to the current, 
see Figs. and d. Secondly, by increasing the bias volt- 
age, one increases the number of nanotube charge states 
participating in transport. This leads to a number of 
peaks in the differential conductance, which can be seen 
in Figs. HJj and e. On the other hand, the corresponding 
Fano factor is displayed in Figs. ^ and f. As already 
mentioned above, the Fano factor diverges in the limit 
of zero bias due to thermal noise. When \eV\ « kBT, 
the shot noise becomes dominant, and the Fano factor 
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FIG. 4; Bias dependence of the current (a,d), differential 
conductance (b,e) and the Fano factor (c,f) for Vg = —0.2666 
V (left panel) and for Vg = -0.1202 V (right panel). The 
parameters are the same as in Fig. [H 



is then significantly reduced. However, it is still super- 
Poissonian in the regions, where the current is exponen- 
tially suppressed (blocked either by the charging energy 
or by the absence of discrete levels available for tunneling 
in the tunneling window). The factor F becomes signifi- 
cantly reduced below the Schottky value, F < 1, for the 
bias voltages large enough to set the system beyond the 
blockade regions. The corresponding Fano factor is then 
slightly above 0.5, which clearly shows suppression of the 
shot noise in comparison to the Poissonian transport pro- 
cesses, and also indicates the important role of Coulomb 
correlations in transport. 

The super-Poissonian shot noise in the Coulomb block- 
ade regime is particularly pronounced close to the thresh- 
old voltage above which the sequential transport is al- 
lowed, see Fig. HI; and f. This enhancement can be 
accounted for as follows. Assume that the voltage is 
slightly below the threshold one (close to the maximum 
of the peak in the Fano factor). The system is then in the 
blockade regime, where the current is exponentially sup- 
pressed. A nonzero small current can flow due to thermal 
excitations. There is an exponentially small probability 
that one electron leaves the CNT (reducing its charge 
state) and then another electron tunnels to the CNT, 
cither to the same energy level or to a level of higher 
energy. If it tunnels to the level of higher energy, trans- 
port through this level is allowed and the electron can 
easily leave the nanotube, while another one can jump 
to the same level or to the level of lower energy. If it 
tunnels to the same level, further tunneling processes are 
allowed. If it tunnels to the low energy level, the system 
becomes blocked again. Similar scenario also holds when 
an electron tunnels first to the CNT (increasing its charge 
state) and then leaves the CNT. All this leads to large 



fluctuations in the current, and consequently to super- 
Poissonian shot noise. However, as we have already men- 
tioned above, the sequential transport is negligible in the 
blockade regions and the current is dominated there by 
higher-order (e.g. cotunneling) contributions (not con- 
sidered here). Thus, the Fano factor in the blockade 
regions may be significantly modified (reduced) by the 
cotunneling current. 




V (mV) 

FIG. 5: Tlie Fano factor, calculated using the total noise 
(solid line) and the excess noise (dashed line), as a function 
of the bias voltage for Vg — —0.1202 V. The parameters are 
the same as in Fig. [1] 

We also notice that experimentally one usually mea- 
sures the excess noise, S* = S{V) — S{V = 0). The corre- 
sponding Fano factor is then defined as, F* = S* / (2e|/|). 
In Fig. O we show the bias voltage dependence of the 
Fano factor F for Vg = —0.1202 V, calculated using the 
total shot noise, S, and F* calculated using the excess 
noise, S*. One can see that in the sequential tunnel- 
ing regime the two definitions yield equal Fano factors, 
F F*. This is due to the fact that the noise due to 
sequential tunneling in the Coulomb blockade regime is 
negligible as compared to the noise above the threshold 
voltage, S{V = 0) <C S{V). Interestingly, the super- 
Poissonian behavior observed in the Coulomb blockade 
is also present irrespective of the definition of the Fano 
factor. The only difference between F and F* can be 
seen in the low bias voltage regime where F oo and 
i^* — > 0, as y — > 0, see Fig. [5] Therefore, in this paper 
we only show and discuss the Fano factor calculated from 
the total shot noise, S, while we notice that the behav- 
ior of F exactly reflects the behavior of F*, except for 
V -^0. 



IV. CARBON NANOTUBES COUPLED TO 
FERROMAGNETIC LEADS 

Carbon nanotubes connected to ferromagnetic 
leads are of particular importance for modern 
magnetoelectronics<^2i^ It is worth noting that the 
magneto-electronics based on all-metal multilayers has 
been making impressive progress over the last two 
decades. Owing to the discovery of the giant and 
tunnel magnetoresistance effects, the multilayer-based 
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FIG. 6: (color online) Differential conductance in the parallel 
(a) and antiparallel (b) magnetic configurations and tunnel 
magnetoresistance (c) as a function of bias and gate voltages. 
The parameters are the same as in Fig. [T] except that pl = 
PR = 0.4. 



devices are now wide-spread in the modern electron- 
ics (nonvolatile memory elements, angular meters, 
read/write heads, etc.). Intensive studies have also been 
carried out on all-semiconductor magnetic systems, 
since the discovery of ferromagnetic semiconductors4i 
The hybrid systems consisting of ferromagnetic metals 
or semiconductors and molecules^i^ open a new area 
in modern magnetoelectronics, and therefore draw 
attention of physicists, engineers and technologists. A 
very promising example of magnetic metal/ molecule 
junctions are ferromagnetically contacted CNTs, studied 
in this section. 

Wc consider thus electronic transport in the situation 
when the CNT is weakly coupled to two ferromagnetic 
leads. The charge current depends then on magnetic con- 
figuration of the system, i.e. on the relative alignment of 
the leads' magnetic moments. Differential conductance 
in the parallel and antiparallel magnetic configurations 
as a function of the bias and gate voltages is shown in 
Fig. [6^ and b. As one can note from the corresponding 
scales, the conductance is smaller in the antiparallel con- 
figuration. This is a general behavior, characteristic of 
the so-called normal spin valves (normal spin valve ef- 
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FIG. 7: Tunnel magnetoresistance as a function of gate volt- 
age for several values of the bias voltage. The parameters 
are the same as in Fig. [T] with pl = Pr = 0.4. Tunnel 
magnetoresistance in the linear response is given by TMR = 
PlPr/(1 -pLPn). 

feet), and results from the spin asymmetry in tunneling 
processes. Apart from this difference, the general behav- 
ior of differential conductance with transport and gate 
voltages in both configurations is qualitatively similar, 
with characteristic diamonds corresponding to blockade 
regions and four-fold periodicity in the linear response 
regime (similar to those found in CNTs contacted to non- 
magnetic leads, see Fig. [5]). 

A quantity which takes into account the difference be- 
tween Fig. [6^ and b, and which is of practical impor- 
tance, is the so-called tunnel magnetoresistance (TMR). 
The TMR phenomenon consists in a change of the tun- 
neling current when magnetic configuration varies from 
antiparallel to parallel alignment, and is quantitatively 
defined 

TMR = ~ , (8) 
Iap 

where Ip (Iap) is the current flowing through the bi- 
ased system in the parallel (antiparallel) magnetic con- 
figuration. The dependence of TMR on the transport 
and gate voltages is shown in Fig. [5J;. From this plot 



one can conclude that TMR in the strictly linear re- 
sponse regime {V — > 0) is constant, i.e. independent 
of the gate voltage. More precisely, it is equal to a half 
of the JuUiere value4^ for the corresponding planar junc- 
tion, i.e. TMR = plPr/(1 — PlPh)- This behavior is 
rather general for electronic transport through quantum 
dots when only first-order transport processes arc taken 
into account, and holds no longer when higher-order pro- 
cesses are included - especially in the blockade regions<^ 
We note again, that including higher-order contributions 
may significantly change the TMR in the blockade re- 
gions, i.e. in the areas of Fig. [St corresponding to the 
black diamonds in Fig. [6^ and b. Generally, one may 
expect some enhancement of TMR in these regions, sim- 
ilarly to the enhancement of TMR due to cotunncling 
processes, observed in other magnetic tunnel junctions 
including also granular structures. 

As follows from Fig. \Et, TMR is no longer constant in 
the nonlinear transport regime. This is shown explicitly 
in Fig. [3 where the gate voltage dependence of TMR is 
presented for several values of the bias voltage. The TMR 
is always positive for the assumed values of the parame- 
ters, and oscillates with the gate voltage (contrary to the 
linear response regime, where TMR is constant). The 
amplitude of the oscillations, however, decreases with in- 
creasing the bias voltage. 

The gate voltage dependence of the Fano factor in the 
parallel and antiparallel magnetic configurations is shown 
in Fig. [5] for a few values of the bias voltage. The gen- 
eral features of the Fano factor with increasing bias and 
gate voltages are similar to those observed in the case of 
CNTs contacted to nonmagnetic electrodes. As before, 
the noise is super-Poissonian in the blockade regions, and 
sub-Poissonian outside the blockade regions. However, 
the Fano factor F in the antiparallel configuration is dif- 
ferent from that in the parallel one, and may be either 
larger or smaller than the latter. This difference is par- 
ticularly pronounced for larger values of the bias voltage, 
where the system is out of the blockade regions. We note 
that for the value of the spin polarization of the leads 
assumed in calculations, the Fano factor in the antipar- 
allel configuration, Fap, is larger than the Fano factor in 
the parallel configuration, Fp . The ratio of Fp /Fap how- 
ever depends on the spin polarization of the leads, and 
for half-metallic electrodes Fp becomes super-Poissonian, 
while Fap approaches unity. This behavior is similar to 
that observed in quantum dotsj ^^i^^'^" 

The bias voltage dependence of the Fano factor is 
shown explicitly in Fig. [5] for two different values of the 
gate voltage chosen so that the first (second) value cor- 
responds to the middle of the small (large) diamond in 
Fig. [51 Figure [H] clearly shows that the shot noise is sub- 
Poissonian, except for the blockade regions, where the 
noise is enhanced above the Schottky value. This be- 
havior is similar to that observed in CNTs contacted to 
nonmagnetic leads, and its physical origin is the same as 
it does not result from magnetism of the electrodes. The 
effect is only quantitatively modified by ferromagnetism 
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FIG. 8: (color online) The Fano factor in the parallel (solid) 
and antiparallel (dashed) magnetic configurations as a func- 
tion of gate voltage for several values of the bias voltage. The 
parameters are the same as in Fig. [T]with pl — Pk = 0.4. 
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FIG. 9: (color online) Bias dependence of the current (a,e), 
differential conductance (b,f), TMR (c,g), and Fano factor 
(d,h) for Vg = -0.2666 V (left column) and for Vg = -0.1202 
V (right column), and Pl = Pr = 0.4. The solid (dashed) line 
in (a,e), (b,f) and (d,h) corresponds to the parallel (antipar- 
allel) magnetic configuration. The other parameters are the 
same as in Fig. [1] 



of the electrodes - the magnitude of the super-Poissonian 
shot noise is different in parallel and antiparallel configu- 
rations, as shown in Figs. [8] and [9l In addition, in Figs. [9^ 
and e we also show the current flowing through the sys- 
tem in both magnetic configurations. The difference in 
currents, resulting from spin asymmetry in tunneling pro- 
cesses in the antiparallel configuration, leads to nonzero 
TMR displayed in Figs.jS}: and g. On the other hand, the 
bias dependence of the differential conductance displays 
characteristic peaks, indicating the different charge states 
of the nanotube being active in transport, see Figs, 
and f . 



V. SUMMARY AND CONCLUSIONS 

Using the real-time diagrammatic approach we have 
analyzed transport through CNTs contacted to nonmag- 
netic and also ferromagnetic electrodes in the weak cou- 
pling regime, where the effects due to charging of CNT 
with individual electrons are clearly seen in transport 
characteristics. Numerical results on the conductance 
in linear response regime reproduce quite well the ex- 



perimentally observed four-peak patterns. Generally, 
conductance in both the linear and nonlinear transport 
regimes reveals features of the discrete electronic struc- 
ture of the CNT-based quantum dots. When the CNT is 
contacted to ferromagnetic leads, the difference in con- 
ductance in the parallel and antiparallel configurations 
leads to TMR, which also reveals characteristic features 
of the discrete charging with single electrons and discrete 
intrinsic electronic structure of the CNTs. 

We have also found that in the sequential tunneling 
regime the shot noise of the single wall metallic carbon 
nanotubes is generally suppressed below the Schottky 
value, except for the blockade regions, where it is larger 
than the Schottky value. For voltages above the thresh- 
old for sequential tunneling the corresponding Fano fac- 
tor has been found to be slightly above 0.5, while in the 
blockade regions it is larger than 1. However, we ex- 
pect that the shot noise in the blockade regions will be 
strongly modified by higher-order contributions to the 
current. It has also been shown that for the assumed 
parameters the Fano factor in the antiparallel configura- 
tion is typically larger than the Fano factor in the parallel 
configuration. 



9 



Acknowledgments Science and Higher Education as a research project in 

years 2006-2009. I.W. also acknowledges support from 
This work was supported by: the EU grant CARDEQ t^ie Foundation for Polish Science, 
under contract IST-02 1285-2, and the Pohsh Ministry of 



* Electronic address: |weymann@amu.edu.pl| 

^ L. P. Kouwenhoven, D. G. Austing and S. Tarucha, Rep. 

Prog. Phys. 64, 701 (2001). 
^ P. Harrison, Quantum Wells, Wires, and Dots: Theoretical 

and Computational Physics (Wiley, Chichester, 1999). 
^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 

Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka, 

and D. M. Treger, Science 294, 1488 (2001). 

* Semiconductor Spintronics and Quantum Computation, 
ed. by D. D. Awschalom, D. Loss, and N. Samarth 
(Springer, Berlin 2002). 

^ S. Maekawa and T. Shinjo, Spin Dependent Transport in 
Magnetic Nano structures (Taylor & Francis 2002). 

® I. Zutic, J. Fabian, S. Das Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

S. Sanvito and A. R. Rocha, J. Comput. Theor. Nanos. 3, 
624 (2006). 

® S. Maekawa, Concepts in Spin Electronics (Oxford 2006). 
^ R. Saito, M. S. Dresselhaus, and G. Dresselhaus, Physical 

Properties of Carbon Nanotubes (London, UK: Imperial 

College Press, 1998). 
^° M. P. Anantram and F. Leonard, Rep. Prog. Phys 69, 507 

(2006). 

W. Hoenlein, G. S. Duesberg, A. P. Graham, F. Kreupl, 
M. Liebau, W. Pamler, R. Seidel, and E. Unger, Micro- 
electron. Eng. 83, 619 (2006). 

K. Tsukagoshi, B. W. Alphenaar, and H. Ago, Nature 401, 
572 (1999). 

" B. Zhao, I. Monch, H. Vinzelberg, T. Miihl, and C. 
M. Schneider, Appl. Phys. Lett. 80, 3144 (2002); J. Appl. 
Phys. 91, 7026 (2002). 

Jae-Ryoung Kim, Hye Mi So, Ju-Jin Kim, Jinhee Kim, 
Phys. Rev. B 66, 233401 (2002). 

S. Sahoo, T. Kontos, J. Purer, C. Hoffmann, M. Graber, 

A. Cottet, and C. Schonenberger, Nature Physics 1, 102 
(2005). 

^® A. Jensen, J. R. Hauptmann, J. Nygard, and P. E. Linde- 
lof, Phys. Rev. B 72, 035419 (2005). 

H. T. Man, I. J. W. Wever, and A. F. Morpurgo, Phys. 
Rev. B 73, 241401 (R) (2006) 

L. W. Liu, J. H. Fang, L. Lu, H. F. Yang, A. Z. Jin, and 
C. Z. Gu Phys. Rev. B 74, 245429 (2006). 

B. Nagabhirava, T. Bansal, G. U. Sumanasekera, and B. 
W. Alphenaara, L. Liu, Appl. Phys. Lett. 88, 023503 
(2006) 

S. Krompiewski, phys. stat. sol. (b) 242, 226 (2005). 

S. Krompiewski, Semicond. Sci. Technol. 21, S96 (2006). 
^2 A. Cottet and M-S. Choi, Phys. Rev. B 74, 235316 (2006). 
23 A. Cottet, T. Kontos, S. Sahoo, H. T. Man, M.-S. Choi, W. 

Belzig, C. Bruder, A. F. Morpurgo, and C. Schnenberger, 

Semicond. Sci. Technol. 21, S78 (2006). 
2'* C. Schonenberger, Semicond. Sci. Technol. 21, SI (2006). 
2^ W. Liang, M. Bockrath, and H. Park, Phys. Rev. Lett. 88, 

126801 (2002). 



S. Sapmaz, P. Jarillo-Herrero, J. Kong, C. Dekker, L. P. 
Kouwenhoven, and H. S. J. van der Zant, Phys. Rev. B 71, 
153402 (2005). 

2^ Y. Oreg, K. Byczuk, and B. I. Halperin, Phys. Rev. Lett. 
85, 365 (2000)'. 

2* R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 

(1997); Eur. Phys. J. B 3, 281 (1998). 
2^ C. Kane, L. Balents and M. P. A. Fisher, Phys. Rev. Lett. 

79, 5086 (1997). 
3° M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. E. 

Smalley, L. Balents and P. L. McEuen, Nature 397, 598 

(1999). 

L. Mayrhofer and M. Grifoni, Phys. Rev. B 74, 121403 

(2006) . 

L. Mayrhofer and M. Grifoni, Eur. Phys. J. B 56, 107 

(2007) . 

33 S. Koller, L . Mayrhofer, M. Grifoni, W. Wetzels, 

jcond-mat/070 3756 (unpublished). 
3'' E. Onac, F. Balestro, B. Trauzettel, C. F. J. Lodewijk, and 

L. P. Kouwenhoven, Phys. Rev. Lett. 96, 026803 (2006). 
35 F. Wu, P. Queipo, T. Tsuneta, T. H. Wang, E. Kauppinen, 

and P. J. Hakonen, cond-mat/0702332 (unpublished). 
3** T. Tsuneta, P. Virtanen, F. Wu, T. Wang, T. T. Heikkila, 

and P. J. Hakonen, cond-mat/0706.2658 (unpublished). 
3^^ D. V. Averin and A. A. Odintsov, Phys. Lett. A 140, 251 

(1989); D. V. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 

65, 2446 (1990); K. Kang and B. L Min, Phys. Rev. B 55, 

15412 (1997). 

3* L Weymann and J. Barnas, Phys. Rev. B 75, 155308 
(2007). 

3^ L Weymann, J. Konig, J. Martinek, J. Barnas, and G. 
Schon, Phys. Rev. B 72, 115334 (2005). 
H. Schoeller and G. Schon, Phys. Rev. B 50, 18436 (1994); 
J. Konig, J. Schmid, H. Schoeller, and G. Schon, Phys. 
Rev. B 54, 16820 (1996). 

A. Thielmann, M. H. Hettler, J. Konig, and G. Schon, 
Phys. Rev. B 68, 115105 (2003); Phys. Rev. B 71 045341 
(2005). 

Ya. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 
■*3 A. Nauen, F. Hohls, N. Maire, K. Pierz, and R. J. Hang , 
Phys. Rev. B 70, 033305 (2004). 

H. Ohno, D. Chiba, F. Matsukura, T. Omiya, E. Abe, T. 

Dietl, Y. Ohno, and K. Ohtani, Nature 408, 946 (2000). 

W. J. M. Naber, S. Faez and W. G. van der Wiel, 

cond-mat / 0703455 (unpublished) . 

M. Julhere, Phys. Lett. A 54, 225 (1975). 

J. Barnas and A. Pert, Phys. Rev. Lett. 80, 1058 (1998). 

B. R. Bulka, Phys. Rev. B 62, 1186 (2000). 

A. Cottet, W. Belzig and C. Bruder, Phys. Rev. B 70, 
115315 (2004). 

5° L Weymann and J. Barnas, J. Phys.: Condens. Matter 19, 
096208 (2007). 



